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qq , Abstract. We briefly review a matrix based method to compute the Casimir operators of 

■ Lie algebras, mainly certain type of contractions of simple Lie algebras. The versatility of 
the method is illustrated by constructing matrices whose characteristic polynomials provide 
the invariants of the kinematical algebras in (3+l)-dimensions. Moreover it is shown, also 

, for kinematical algebras, how some reductions on these matrices are useful for determining 

■ the missing operators in the missing label problem (MLP). 
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1 Introduction 



Casimir operators of Lie algebras were originally introduced in the frame of representation 
theory of semisimple algebras, and were soon recognized as a powerful tool for the structural 
analysis. They were determined by means of the universal enveloping algebra of a Lie algebra, 
and successive applications of Lie groups and algebras to other mathematical and physical 
problems led to alternative methods to compute these invariants, like the orbit method or the 
formulation in terms of differential equations, which allow to generalize the concept of invariant 
functions on Lie algebras beyond the pure algebraic frame. 

In 1950 I.M. Gel'fand showed how to use the generic matrix of standard representation of 
orthogonal groups to determine the Casimir operators of the algebra by means of characteristic 
polynomials. More specifically, for the pseudo-orthogonal Lie algebra so(p, q) (N = p + q) given 
by the ^N(N — 1) operators E^ u = — E vfl and having brackets: 

[E/ju,, Exu] = g^\E ua + g^Exu — gy\E^ a - g^Ex^, 

where g = diag (1, . . . , 1, —1, . . . , —1), the Casimir operators were obtained with the formula pQ: 

N 

N-k 



P(A) = \M M - \U N \ = X N + C k \ 



k=l 
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M Pt g being the matrix 
/ • 



M, 



p.q 



~9jj e ij 







-Qnnzin 



-QNN&jN 



\ e\N ■ ■ ■ 9jjdjN ■■■ / 

and Id^r the identity matrix of order N. Once symmetrized, the functions Ck provide the classical 
operators in the enveloping algebra of so(p,q). This direct matrix approach seems natural, 
taking into account that the eigenvalues of Casimir operators are essential for the labelling of 
representations of the Lie algebra. By application of the classical theory of semisimple Lie 
algebras, similar formulae can be developed for other semisimple Lie algebras j2J. It is natural 
to ask whether the essence of the Gel'fand method can be generalized to other non-semisimple 
Lie algebras, by means of extended matrices that possibly correspond to some representation of 
the algebra. The two main interesting cases are: 

1. Generalized Inonii-Wigner contractions of (semi)simple Lie algebras. We have inhomoge- 
neous algebras as a special case. 

2. Semidirect products of simple and Heisenberg Lie algebras. 

The method usually employed to determine the invariants of a Lie algebra consists in solving 
of a system of partial differential equations (PDEs) related to a realization of the Lie algebra 
as differential operators [HI 13- If {Xi, • • • ,X n } is a basis of q and {C^} is the structure tensor, 
the realization of q in the space C°°(g*) is given by: 



Xi 



-CijXkd Xj , 



where [Xi, Xj] = C^X^ (1 < i < j < n, 1 < k < n), {x\, . . . , x n } is the corresponding dual basis, 
q* being the dual space to q. Invariants of the coadjoint representation are functions on the 
generators F (X±, . . . , X n ) of q such that [Xi,F (X\, . . . , X n )} = 0, and are determined by solving 
the system of PDEs: 



X l F(x l , 



j Xn J 



dF 



■CjjXfc — (x\, . . . , x n ) — 0, 



1 < i < n. 



(2) 



Classical Casimir operators are recovered with replacing the variables x% by the corresponding 
generator Xi (possibly after symmetrizing). The cardinal M (q) of a maximal set of independent 
solutions is described in terms of the following formula: 

M(g)= dim - rank A{q), (3) 

where A (q) is the matrix representing the commutator table of g over a given basis, i.e., 

A( fl ) = (Cgx fc ). 

In this work we briefly review some recent work on matrix procedures that generalize the 
Gel'fand formula for various types of Lie algebras, specifically inhomogeneous Lie algebras. 
Applications to the computation of missing label operators are also given. 



2 Contractions of simple Lie algebras 

Let q be a simple Lie algebra of rank k such that the Casimir operators of q are given by the 
characteristic polynomial P (T) of a matrix 

A = Y,Fi {Xi) Xi , 
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where Ti is equivalent to the iV-dimensional standard representation of 0. Suppose further that 
the Lie algebra q' is a contraction of q and that M(g) = M(g'). To compute the Casimir operators 
of g' using a matrix extension of A, the following procedure by steps has been proposed 0: 

1. Take the matrix A = Ti (JQ) Xj. 

2. Let \t € be the automorphism giving the contraction 1 and consider the matrix A' = 
X^jTx {X'i)) x'i giving the invariants of g over the transformed basis \X[ := $ e (Xj)}. 

3. Let P e (T) := \ A' -TU N \ and a := deg e P £ (T). Take the limit P(T) := lim iP e (T). 

4. Rewrite -P(T) as a linear combination of determinants. 

5. Define Ij = Cj as the polynomial coefficients of P (T). 

6. Check the functional independence of the functions Ij. 

7. If necessary, symmetrize the functions Ij in order to recover the Casimir operators of g' . 

The validity of the method is formally proved using the contraction explicitly. If ^ e is the 
automorphism of defining the contraction, then the matrix 

A' k = J2n (xi) x\ 

i 

expresses the Casimir operators of over the transformed basis. Developing the corresponding 
characteristic polynomial P(T) and taking into account the limit, after some algebraic manip- 
ulation we can obtain an expression of P(T) that does not involve the contraction parameter 
e anymore. By transitivity of contractions, the same algorithm can be applied formally to the 
contraction of reductive algebras. 

For Lie algebras of the type ws = s~®rhN, i-e., semidirect products of a simple and a Heisen- 
berg Lie algebra, the argument remains valid for the contraction ws s0r(2iV + 1)L±. The 
structure of representations of simple Lie algebras compatible with Heisenberg algebras was an- 
alyzed in Depending on this structure, we obtain by contraction either an inhomogeneous 
Lie algebra or a general affine algebra. The first step is to compute invariants of ws. This is 
done by constructing a copy of the Levi part 5 in the enveloping algebra of ws 2 , which can be 
used as the simple algebra to which the preceding procedure is applied. 

Remark 1. In general, the shape of P(T) will depend essentially on the structure of the standard 
representation and the contraction considered, and possibly involves matrices depending on T. 
In particular, it will be different for the cases where this representation is by real or complex 
matrices, i.e., if it is of first or second genus [H]. This fact can also originate some dependence 
problems. 

2.1 Inhomogeneous symplectic algebras 

As an example of a general class to which the algorithm applies, we consider the inhomogeneous 
symplectic Lie algebras Isp(2N, R) and the matrix formula obtained for them in Over the 
basis {Xi t j,X-ij,Xij,Pi,Qi} (1 < i, j < N) the brackets of Isp (2N,M) are given by: 

[Xij,Xk,i] = SjkXu — 5uXkj + EiEjdj-iXk-i — EiEjdi-kX-j^, 
[Xij,Pk] = SjkPi, [Xij,Qk] = —SikQj, 

Qk] = —SjkPi — $ikPj, [Xi,-j, Pk] = $ikQj + $kjQi- 

1 By this we mean that the limit [X, Y] := lim [* e (X), * e (K)] exists for any X, Y G g. This equation 

e — >oo 

defines a Lie algebra @' called the contraction of g (by ^ e ). 

2 These variables can be found analyzing the noncentral Casimir operator of the subalgebras s'©rfj(iV), where 
s' is a simple subalgebra of rank one of s 
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Table 1. Nonisomorphic kinematical algebras in (3 + 1) dimensions |1U| . 





so (4,1) 


so (3, 2) 


Iso (3,1) 


Iso (4) 


Ne exp 


Nc osc 


Carroll 


G(2) 


Static 


[B,P] 


K 


K 





K 


K 


K 











[H,K] 


P 


P 


P 





P 


P 





P 





[P,P] 


J 


-J 





J 

















[K,K] 


J 


-J 


J 




















[P,K] 


H 


H 


H 


H 








H 









Proposition 1. Let N > 2. Then the Casimir operators C 2 k of Isp (2N,~. 
coefficients of the polynomial 



are given by the 



N 



\C - TId2AT+i| + |C2JV+l,2iV+l - 



2k 



where 



( x h 



C 



XN,l 
Xl-1 



xn,n —x-i 7 n 

X\-N -XlA 



Xl-N ■ ■ ■ XN-N —Xl,N 

\ -qi ■■■ -q N pi 



k=l 

-X-l >N piT \ 

~X-N,N PnT 

xn,i qiT 

—xn,n <1nT 

PN / 



and C2N+i,2N+i is the minor ofC obtained deleting the last row and column. Moreover deg C 2 k+i 
2k + 1. 



3 The kinematical algebras in (3 + l)-dimensions 

As an interesting physical example, we apply the preceding procedure to the class of kinematical 
Lie algebras in (3 + l)-dimensions. This choice is appropriate since it contains algebras of both 
types and because they are all related by contractions thus by transitivity the procedure 

can be applied. Although their invariants have been obtained repeatedly in different contexts |12| 
I13j . it is worthy to be done on the basis of the above arguments, which, moreover, show that 
the matrix providing the invariants is not necessarily related to a faithful representation of the 
algebra. 

Following the original notation of JO]) kinematical Lie algebras are defined over the basis 
i<«<3' where J a are spatial rotations, P a spatial translations, K a the boosts 
and H the time translation, constrained to the condition of space isotropy 

[Ja,M =e Q/3 V 7 , [J a ,P^]=e Q ^P 7 , [J a ,Kp] = e a ^K 1 , [J a ,H}=0, 

as well as the assumption that time-reversal and parity are automorphisms of the group. Taking 
the compact notation [X,Y] = Z for [X a ,Yg] = E a @' y Z~ / , the brackets of the nonisomorphic 
kinematical Lie algebras are given in Table l 3 . 

As shown in ^01 > any kinematical Lie algebra is obtained by contraction of the simple de Sitter 
Lie algebras so (3, 2) and so (4, 1). With the exception of the static algebra, all of the remaining 
ones possess two independent invariants. 

3 We have omitted the Para-Poincare and Para-Galilei Lie algebras, since they are isomorphic to the Poincare 
and Galilei algebras, respectively, although they are physically different. For purposes of invariants, this physical 
distinction is irrelevant. 
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3.1 De Sitter algebras 



Since both de Sitter algebras are simple and pseudo-orthogonal, their Casimir operators follow at 
once from application of the Gel'fand formula. For obtaining the invariants over the kinematical 
basis above, the matrix Q has to be slightly transformed. 

1. Anti de Sitter algebra so (3, 2) . The matrix related to the standard representations is: 



D 



( J3 

-h o 

-h -ji 

—ki fa 



32 

ji 

-fa 



-fa pi \ 

fa ~P2 

-fa P3 

h 

-h / 



(4) 



\ Pi ~P2 P3 

Computing the characteristic polynomial we have \D — TTdsl = T 5 + C2T 3 + C4T, where 
C 2 



jaj a ~ PaP° ~ fak a + h 
,-ai.2 



Ci = j a j a h 2 + {( Pa p a )(k a k a ) - ( Pa k a ) 2 ) - (j a p a f - (j a k a ) 2 - 2e a ^j a p^h. 
2. De Sitter algebra so (4, 1). The resulting matrix is similar to the previous one: 



D 



( J 3 

-h 

-32 -h 

-fa k 2 

V -Pi P2 



h -fa Pi \ 

ji fa -P2 

-fa P3 

-fa h 

-P3 h / 



Then we have \D - Tld 5 \ = T 5 + C 2 T 3 + C 4 T, where 
C 2 = j a j a + PaP a -k a k a -h 2 , 

Ci = -j a j a h 2 - {( PaP a )(k a k a ) - ( Pa k a ) 2 ) + (j aP a ) 2 - (j a k a ) 2 + 2e al3 y aP ^fah. 
For these two algebras, the result follows at once from the Gel'fand formula. 



3.2 The nonrelativistic cosmological Lie algebras 

The Newton algebras Ne + and Ne~ are obtained as contractions of the de Sitter and Anti de 
Sitter algebras, respectively. Since the subalgebra generated by {K a ,P a } is Abelian, it will 
follow that the invariants of these algebras will not depend on the rotation and time-translation 
generators. 

1. The Newton algebra Ne~ [oscillating universe]. It can be easily verified that this Lie 
algebra is indeed an extension by a derivation of the nine dimensional Lie algebra 
so (3) ©2a(teo(3)6£i- By the comment above, the matrix giving the Casimir operators of 
Ne~ will not be related to a faithful representation of the algebra. We have: 



/ 











-fa 


Pi 


\ 













fa 


~P2 















-fa 


P3 






-fa 


fa 


-fa 










V 


Pi 


~P2 


P3 








/ 



Expanding the secular equation we arrive at \D — TId§\ = T 5 + C2T 3 + C4T, where 
C 2 = -p a p a -k a k a , 

Cl = ((PaP a )(fak a ) - ( Pa k a ) 2 ) ■ 
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2. The Newton algebra Ne + [expanding universe]. This Lie algebra is also an extension by 
a derivation of so (3) ©2adso(3)^i- in this case the matrix to be used is: 



/ 











-h 


Pi 


\ 













k 2 


-Pi 















-ka 


P3 






-ki 


k 2 


-h 










V 


-Pi 


P2 


~P3 








/ 



Then we have \D - Tld 5 \ = T 5 + C 2 T 3 + C 4 T, where 

C 2 = p a p Q -k a k a , 

Ca= {(PaP a )(kak a ) - ( Pa k a ) 2 ) . 



3.3 The inhomogeneous (pseudo)-orthogonal algebras 

In order to obtain the matrix for the inhomogeneous algebras Iso (3, 1) and Iso (4), we use the 
contractions so (3, 2) iso (3, 1) and so (4, 1) ~> Iso (4), respectively. 

1. The Poincare Lie algebra Iso (3, 1). The matrix D is given by 



D 



( 

-J3 

-32 

-h 
V pi 



33 


-31 
k 2 

~P2 



,72 

h 


-fc 3 

P3 



-*!3 



-/t 



Pi \ 

P3 
h 

o / 



We define P (T) = \D — TIds| + T — TM4I, where -D55 is the minor of D obtained 
by deleting the fifth column and row. In particular, the matrix D55 corresponds to that 
of the Lorentz algebra so (3,1). Expanding the expression for P(T), we get P (T) = 
T 5 + C 2 T 3 + C 4 T, where 



C 2 = h z - PaP c 



C4 = j a j a h 2 + (p aP a ) (k a k a ) - (j aP a Y - (p a k c 
Moreover, the matrix D can be decomposed as 



D = D 1 +D 2 



2e a ^ jaPpk^h. 



( 





33 


J2 




Pi 


\ 
















\ 




-33 





ii 


^2 


-P2 























-32 


-31 





-*3 


P3 




+ 



















-h 


k 2 


-h 





h 





















\ 





















\ Pi 


-P2 


P3 


-h 


J 



where D\ defines a faithful representation of Iso (3, 1). 

2. The inhomogeneous algebra Iso (4). Here the polynomial is P (T) = \D — TIds| + T\D^ — 
TM4I, where -D55 is the minor of D obtained deleting the fifth column and row. The 
matrix D is given by 



D 





-33 
-32 

-h 
-pi 



33 


-31 

k 2 

P2 



32 

31 


-h 

~P3 



-h 

k 2 
-k 3 



/; 



Pi \ 

~P2 
P3 
h 

/ 
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Expanding P (T), we get P (T) = T 5 + C 2 T 3 + C 4 T, with 



C 2 
C 4 



j a fh 2 + (p Q p Q ) (A^ Q ) + (j a k a f 
In this case, D decomposes as 



(p a k a ) 2 - 2e a ^j aPp k 1 h. 



( 





33 


32 





Pi 


\ 




f 








-h 


o \ 




-33 





31 





~P2 















k 2 







-32 


-31 








P3 




+ 











-k 3 







-h 


k 2 


-** 





h 





















\ 


-Pi 


P2 


-P3 




















h 


o / 



D = D l +D 2 



and Di is the matrix related to a the faithful representation of Iso (4) by 5 x 5 matrices. 
3.4 The Carroll Lie algebra 

Among the classical kinematical Lie algebras, the Carroll Lie algebra is the only isomorphic to 
the semidirect product of a simple Lie algebra (the compact algebra so (3)) and a Heisenberg 
Lie algebra. Indeed the noncentral Casimir operator can be determined using the determinant 
procedure developed in [S]. However, the Casimir operators (the second is trivial, since the 
centre is nonzero) can also be obtained by the same method as before. 
We define P (T) = \D — TIds|, where D is the matrix 



D 



Observe that in this case, the matrix D is dependent on the variable T. Expanding, we get 
P (T) = T 5 + C 2 T 3 + C 4 T, where 

C 2 = h 2 , 

Ci = Ja j a h 2 + ( PaP a ) (k a k a ) - ( Pa k a 
The matrix D decomposes in this case as 








33 


32 


-h 


Pi 


\ 




~33 





3i 


k 2 


~P2 






~32 


-3i 





-k 3 


P3 






-h 


k 2 


-k 3 


T 


h 




I 


-Pi 


P2 


~P3 


h 


T 


) 



le^jap^h. 



( 





33 


h 





Pi 

















-h 


o \ 




-33 





h 





~P2 

















k 2 







-32 


-3l 








P3 




+ 













-h 







-ki 


k 2 


-h 





h 

















T 





V 





















I 


-Pi 


P2 


~P3 


h 


T J 



D = D 1 +D 2 



Again, the matrix D\ gives rise to a faithful representation of the Carroll algebra. 
3.5 The Galilei algebra G{2) 

As happened for the Newton algebras, the Casimir operators of the Galilei algebra do not depend 
on the variables {j a , h}. Here we consider the polynomial P (T) = \D — TM5I, where 



D 







-h 
-pi 







k 2 

P2 







-ks 

~P3 



-ki 

k 2 
-k 3 







Pi \ 

~P2 
P3 


T 
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Also in this case, the matrix D is dependent on the variable T. Developing the polynomial we 
obtain P (T) = T 5 + C 2 T 3 + C 4 T, where 

C 2 = p a p a , 

C 4 = (PaP a ) (k a k a ) - ( Pa k a ) 2 . 

By the remark above, the preceding matrix is not related to a faithful representation of the 
algebra. 



3.6 The static Lie algebra 



This algebra is nothing but the splittable affine Lie algebra (50 (3) 
commented, it has four invariants, all of the degree two, 



} 2adso{3) 



6L1) 



L As 



h = h, 



PaP 



h = k a k a , h = k a p a . 



To obtain them in matrix form, we consider 



D 



( 





—hi k 2 

\ -Pi P2 



-h piT \ 

k 2 -p 2 T 

-k 3 p 3 T 

-k 3 -hT 

-Pa -h / 



and obtain P (T) = T 5 + (l 2 - if) T 4 - J 3 T 3 - (l 2 I 3 - if) T 2 . Simplifyng the coefficients we 
recover the basis of invariants above. Since the variables associated to the rotations J a do not 
appear in the invariants, the matrix D does not provide a representation of the static algebra. 
For later use we consider the following functions: 



h = h, I 2 =p a p a , h 
h=j a k a , I 7 =jaP a , M 



-.a/37 



h = jaf 
jaPpky. 



h 



k a p a ) 



(5) 



4 Applications: The missing label problem 

As known, irreducible representations of a semisimple Lie algebra are labelled unambigously 
by the eigenvalues of Casimir operators. In a more general frame, irreducible representations 
of a Lie algebra g are labelled usingby means of the eigenvalues of its generalized Casimir 
invariants |14l I15j . The number of internal labels needed is 

i = -(dimfl-A^(fl)). 

If we use a subalgebra f) label the basis states of g, then we obtain i(dimf) + A/"(h)) + I' labels 
from f), where I' is the number of invariants of g that depend only on variables of the subalge- 
bra f) ^3 • In order to label irreducible representations of g uniquely, it is therefore necessary to 
find 

n= -(dimg-A/"(g)-dimf)-A/"(f)))+Z / (6) 

additional operators, which are usually called missing label operators. They are found by integra- 
ting the equations of system (J2J) corresponding to the subalgebra generators. The total number 
of available operators of this kind is easily shown to be m = 2n. 



Application of the Gel'fand Matrix Method to the Missing Label Problem 



9 



In this situation, it is plausible to think that whenever the Casimir operators of a Lie algebra q 
can be determined using determinants of (polynomial) matrices, the same procedure could hold 
for computing missing label operators according to some subalgebra f). In this section we analyze 
the missing label problem for the chain 

so(3) ^ Q, 

where q is a kinematical Lie algebra in (3 + l)-dimensions and so (3) the compact subalgebra 
generated by the {</^}. The missing label operators are among the solutions of the equations: 



dF dF dF dF OF dF 

h— J2WT- +P3T. P2~ h h— h- 

OJ3 Op 2 



dF OF 

dp 3 dk 2 



0, 



oj2 oj 3 dp 2 dp 3 dk 2 dk 3 

dF dF dF dF dF , dF 

-h^~ + h^— -P3^ Vpi-^. k 3 — + h— = 0, 

oji dj 3 dpi dp3 dk\ dk 3 

dF dF dF dF dF dF 

32— Ji^— +P2T> Pin ^ k 2 —, ki-^r- = 0. 

dj 2 dpi 



dp 2 



•diet 



dh 



Due to the space isotropy condition, the above system is valid for any kinematical Lie algebra. 
According to formula there are 

n= X - (dim - M ( ) - dimso (3) - M (so (3))) + l' = ^ (6 - M(q)) + I' 

missing labels. In any case we have I' = 0. Moreover, for any kinematical algebra, with the 
exception of the static algebra, we obtain n = 2 and four available missing label operators, while 
for the static algebra we get n = 1 and m = 2. By using of the matrix notation, the system can 
be rewritten as: 






33 


~32 





P3 


~P2 





h 


-k 2 





-J3 





h 


~P3 





Pi 


-k 3 





ki 





32 


-3i 





P2 


-Pi 





k 2 


-h 









dF \ 



( 9j a 

dF 

dp a 

dF 

dk a 

dF 
dh 



0. 



(7) 



V 9F I 

\ dh / 



Since the matrix has rank three, there are seven independent solutions of the system. The 
number of solutions that do not depend on the variables {j a } of the subalgebra 50 (3) is given 
by [S]: 



J\f' = 7 — rank 






P3 


~P2 





k 3 


-k 2 





~P3 





Pi 


-k 3 





h 





P2 


-Pi 





k 2 


-h 









4. 



It is straightforward to verify that a complete system of independent solutions is given by: 



{h = h, I 2 = p a p a , I 3 = k a k a , I A = j a j a , I 5 = k a p a , Iq = j a k a , I 7 = j a p a } . 



(8) 



In particular, the invariants I\, I 2 , I 3 , I§, which are the independent solutions not involving the 
variables j a , constitute a set of solutions for the static Lie algebra. The function M = e a ^ jaPpk^ 
of ((SJ) is functionally dependent on the previous functions, as shown by the relation 



M 



4h + I7I3 



I 2 I 3 h ~ 2I 5 hI 7 . 



To see how the matrices used for the Casimir operators of kinematical algebras can also be 
applied to the MLP, we consider the Anti de Sitter algebra so (3, 2). In the notation of ©, the 
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Algebra g 




so (4,1) 



50 (3,2) 



iso (4) 



Iso (3, 1) 




Nc" 



Nc+ 



Carroll 





{h,I 2 ,I 6 ,I 7 } {h} 
{h,I 2 ,I 6 ,I 7 } {I 2 } 



Galilei 



Static 



invariants of the algebra are given by C2 = if — h — h + h and C4 = If I 4 + ^2^3 — 15 — — Ij — 
2I\ M, while I4 clearly represents the Casimir operator of so(3). We now look for those missing 
label operators that depend only on the variables {p a ,k a ,h}. To this extent, we consider the 
matrix @ used to compute C2 and C4 and replace the variables j a by 0. We obtain 



/ 











-h 


Pi 


\ 













k 2 


~P2 















-k 3 


P3 






-h 


h 


-h 





h 




V 


PL 


-P2 


P3 


-h 





/ 



Considering the characteristic polynomial we have 



It can be easily verified that {12,-^3,-^5} are solutions of (J7J) independent of {02,04,14}, while 
if is not an independent solution. Therefore the matrix D' provides three of the four available 
missing label operators. The fourth, which can be chosen as 1%, cannot be obtained using D' , 
since it depends on the rotation generators. 

A similar argument can be used for the remaining kinematical algebras q. We consider the 
matrix giving the invariants of q and replace the subalgebra generators j a by 0. Then we 
compute the corresponding characteristic polynomial and see how many independent solutions 
from {02,04,14} are obtained, where C2 and C4 are the quadratic and fourth order Casimir 
operators of q, respectively. Only for the static Lie algebra this method of generating missing 
label operators fails, since their Casimir operators are a maximal set of solutions of system (J7J) 
not depending on the j a . The corresponding results are presented in Table 2. 

5 Final remarks 

The approach presented here to compute Casimir invariants of Lie algebras tries to extend the 
classical results established for semisimple algebras to their contractions, making use of the 



P (T) = \D' - TId 5 1 = T 5 + (If -I 2 - I 3 )T 3 + (hh ~ I$)T. 
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standard representation of the contracted algebra. Although it provides in many cases closed 
formulae of the invariants of contractions, the application of the Gel'fand formula is certainly 
only of interest when the contraction has the same number of invariants than the original algebra. 
Although in the algebras analyzed here no dependence problems have been encountered, one of 
the unsolved problems is to find sufficiency criteria to ensure that the contraction of independent 
invariants of an algebra provides also independent operators of the contraction. Work in this 
direction is in progress. 

As concerns applications, we have seen that the MLP can be analyzed via the generalization 
of the Gel'fand method, by simply reducing the matrices by zeros corresponding to generators 
of the subalgebra considered. This point of view could also be interesting in combination with 
problems in symmetries of differential equations related to contractions of Lie algebras, such as 
the separation of variables |17j . 
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